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Abstract 

The total energies and various bound state properties of the excited 2 1 S(L = 0)— states in 
two-electron helium atoms, including the °°He, 4 He and 3 He atoms, are determined to very high 
numerical accuracy. The convergence of the results obtained for some electron- nuclear and electron- 
electron expectation values and, in particular, for the electron-nuclear and electron-electron cusp 
values, is discussed. The field component of the isotope shift and lowest order QED correction are 
estimated for the 2 1 S(L = 0)— states in the 4 He and 3 He atoms. We also apply our highly accurate 
methods to numerical computations of the excited n 1 S— states (for n = 3 and 4) in two-electron 
atomic systems. 
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I. INTRODUCTION 



In this communication we consider the bound state properties of the excited 2 1 S(L = 
0)— states in two-electron helium atoms: °°He, 4 He and 3 He. In our previous work |l| we 
have performed highly accurate computations for the ground = 0)— state and for 

the triplet 2 3 S(L = 0) -states in a number of helium-like ions. Another work [2j contains 
highly accurate results for the singlet 2 l P(L = 1)— and triplet 2 3 P(L = 1)— states in helium 
atom(s). In this study we want to consider the excited 2 1 S(L = 0)— states in two-electron 
helium atoms. Our current interest to the excited 2 1 S(L = 0)— states in the He atoms stems 
from the following. First, many bound state properties of the excited 2 1 S(L = 0)— states 

n n 

in helium atoms have not been computed in earlier studies (see, e.g., |3(, |4| and references 
therein) . 

Second, there is a common believe that construction of highly accurate wave functions for 
the excited states in two-electron atoms and ions is a significantly more complicated process 
than in the case of the lowest states with the same L and S quantum numbers. Here and 
below the notations L and S designate the quantum numbers of total angular momentum and 
electron spin, respectively (see, e.g., Q and Q). Most of the expected complications most 
likely arise from various numerical instabilities which become crucial during optimization 
of the non-linear parameters in the wave functions. Indeed, the repeated optimization and 
re-optimization of the non-linear parameters in the wave function is a potentially unstable 
process for the excited states. Formally, by optimizing these parameters we are trying to 
decrease the second, third, . . ., n — th eigenvalue, i.e. the second, third, etc, root of the 
eigenvalue equation. The orthogonality of the wave functions of the excited states to the 
wave function of the ground state is not checked during such an optimization. This means 
that after a few steps of the optimization of non-linear parameters the process of construction 
of a highly accurate wave function for an excited state may begin to converge to the ground 
state wave function. This 'instability' problem explains a relatively modest progress achieved 
in highly accurate calculations of the excited states in comparison to the ground state. 

Third, it is commonly assumed that, for some unexplained reasons, the overall conver- 
gence rate for bound state properties of the excited atomic states will be substantially lower 
than analogous convergence rate for the same properties determined for the lowest energy 
state with the same L and S numbers. This problem has not been studied carefully in earlier 
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works. It can be answered by computing a large number of bound state properties for the 
excited 2 1 S'— state in the He atom(s) with the use of different number of basis functions. 



II. METHOD 

Our computational goal in this study is to determine the highly accurate solutions, i.e. the 
eigenstates and corresponding wave functions, of the non-relativistic Schrodinger equation 
#^(ri,r 2 ,r 3 ) = E^(r 1 ,r 2 ,v 3 ), where E < [7|, for 2 1 S(L = O)-state(s) in the neutral 
helium atom. For an arbitrary two-electron atomic system with the nuclear charge Qe the 
non-relativistic Hamiltonian is written in the following form 

H = — — V? - —Vl - -^-(^)Vg - — - — + — , (1) 

2m e 1 2m e 2 2m e ^ M > 6 r 32 r 31 r 21 

where =\ r j — r, |= Tji (i ^ j = (1, 2, 3)) and is the radius- vector of the i— th particle 
written in Cartesian coordinates. Also, in this equation Vj = fgf:, is the gradient 

operator for the particle with index i (i — 1, 2, 3). Here and everywhere below in this 
study the subscripts 1 and 2 designate the two electrons, while the subscript 3 stands for 
the nucleus. In the Hamiltonian, Eq.flJJ), Q = g 3 is the nuclear charge, while M = m 3 3> 1 is 
the nuclear mass. In atomic units where h = l,e = l,m e = 1 the Hamiltomian H, Eq.flT]), 
takes the form 

*=-^-^-^-£-£ + £- < 2 > 

In this study the unknown wave functions are approximated with the use of exponential 
variational expansion in relative coordinates r 32 ,r 31 and r 21 . For the excited 2 1 S(L = 
0)— states, the exponential variational expansion for the spatial part of an arbitrary two- 



electron wave function takes the form 



,8 



1 



^lm = —7=0- + P2i)^2C i exp(-a i u 1 - f3iu 2 - %u 3 ) (3) 
V2 i=1 

1 N 

(-«i^32 - Ann - 7» r 2i) 



(l + P 21 )J2Qexp 
V2 i=1 



where Cj are the linear (or variational) parameters, ojj, and ji are the non-linear parame- 
ters (i = 1, 2, . . . , N) of variational expansion, Eq.fl3]). The operator P 21 is the permutation 
of the two identical (1 and 2) particles (electrons) in the symmetric two-electron ions. In 



this equation u±, 1*2, ^3 are the three perimetric coordinates which are truly independent and 
simply related to the three relative coordinates: U{ = + — r^). The inverse relation 
takes the form = Ui + Uj, where (i,j,k) = (1,2,3). Note that each of the perimetric 
coordinates varies between and +00. 

In general, by using the variational expansion, Eq.([3]), it is possible to construct extremely 
accurate wave functions for the S(L = 0)— states in arbitrary three-body systems, if (and 
only if) the non-linear parameters Qjj, and ji in Eq.(|3]) are varied. In this study we shall 
use the numerical methods and optimization strategy developed in our earlier work [8] . The 
procedure developed in |8|] allows one to produce extremely accurate (or essentially exact) 
variational results for arbitrary three-body systems, including the Ps~ and H~ ions 
Many of these systems, however, have either one ground bound l 1 ^— state, or a few bound 
states with different L and S quantum numbers. For two-electron atoms and ions which 
have infinite numbers of excited states we modified our strategy developed in [8J to provide 
a better numerical stability during each step of optimization. 



III. PROPERTIES 

The results of our variational calculations for the excited 2 1 S(L = 0)— state in the helium 
atom with infinitely heavy nucleus, i.e. in the °°He atom, can be found in Table I which 
contains the total energies E and (r 2 i), (<5 3 i) and z/31 expectation values expressed in atomic 
units. The notations (r 21 ), (£31) and v 3 \ are explained below. These expectation values have 
been determined with the use of trial wave functions with the different number of basis 
functions N. In this study to approximate the highly accurate wave functions we have used 
the exponential variational expansion in relative coordinates r 32 ,r 3 i and r%i, Eq.(j3]). 

Table II contains the expectation values of many atomic properties of the 2 1 S— state 
(expressed in atomic units or a.u.) determined using the most accurate wave functions 
obtained in this study. The physical meaning of almost all of the expectation values pre- 
sented in Tables I and II is clear, and we make here only a few following remarks. The 
expectation values of interparticle distances and their powers are designated as (rfj) for 
k = 1, 2, 3, 4; — 1, — 2 and (ij) = (31), (21) for the two-electron He atom. The expectation 
values of the electron-nuclear, electron-electron and triple delta-functions are designated as 
(S(r 31 )) = (6 31 ),(6(r 21 )) = (5 21 ), (5(r 32 i)) = (5(r 32 )5(r 3 i)) = (6 321 ). Also, Tables I and II 
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include the electron-nuclear (z/ 3 i) and electron-electron (z/ 2 i) cusp values [10|, [11] : 



v ^^t (4) 

where Sfaj) = is the appropriate delta-function and (ij) = (31), (21). In this study we 



omb potentials. 



10j . Ill ], where 



shall assume that all point particles interact with each other by the Cou 
Therefore, the expected (exact) two-particle cusp equals i^j = qiqj ™*™^. 
qi,qj are the particle's charges and m^mj are their masses. For the considered singlet 
2 1 S(L = 0)— states in the two-electron helium atom with the infinitely heavy nucleus °°He, 
the expected electron-nuclear cusp equals — Q = —2, while the electron-electron cusp equals 
0.5. For the 3 He and He helium atoms the electron-electron cusps also equal 0.5 (exactly), 
while electron-nuclear cusps include some small correction for the ?jf — -h ratio. For the 
3 He and 4 He helium atoms one finds for the electron-nuclear cusp (in atomic units) 

v = -Q-^— = -Q(l —) (5) 

where M ^> 1 and Q is the nuclear electric charge. In our calculations of the 2 1 S— states in 
the 3 He and 4 He helium atoms the followin g va lues of nuclear mass were used: M( 3 He) = 
5495.8852 m e and M( 4 He) = 7294.2996 m e [12|, [13|. 

The coincidence between the expected and computed cusp values provides a very con- 
venient, accurate and universal criterion to control the overall quality of the wave function 
constructed in our computations. In general, the electron-electron cusp is a better criterion 
of the overall quality of the variational wave function than the electron-nuclear cusp. This 
means that in actual calculations it is much harder to obtain a good/excellent agreement 
for the electron-electron cusp, than a similar agreement for the electron-nuclear cusp. 

Table II also contains the expectation values of the two interparticle cosine— functions 
which are determined as follows 

Tij = (cos(r ifc A r ifc )) = ( r ' fc Tjk ) , (6) 

where (i,j,k) = (1,2,3). The expectation values are always < 1. The absolute value of 
T2i(= T12) can be considered as an 'ideal' measure of the static electron-electron correlations 
in the two-electron atomic systems. Let us define the quantity (/) which is expressed in 
terms of the relative coordinates (r 3 i, r 32 , r 2 i) or perimetric coordinates (u\, U2, U3) as follows: 

(/) = 7i(if> I — — I xj}) = III I ^( Ul ,u 2 ,u 3 ) | 2 u 1 u 2 u 3 du 1 du 2 du 3 . (7) 
2 r 32 r 3 i r 2 i J J J 



It can be shown that the equality 



7-21 + r 32 + r 3 i = 1 + 4(/> (8) 

holds for arbitrary three-body system. For the two-electron (i.e. symmetric) ions/atoms we 
always have T32 = T31, and therefore, t%\ + 2t 3 i = 1 + 4(/). It can be also shown that in 
an arbitrary Coulomb three-body system < (/) < 0.085. The (/) value can be calculated 
either directly from Eq.([7j), or by applying the expectation values of the cosine functions re- 
computed earlier. The coincidence of these two values of (/) indicates that the t 3 2,t 31 ,T2i 
and (/) expectation values have been computed correctly. 
The virial factor 77 in Table II is determined as follows: 

where (T) and (V) are the expectation values of the kinetic and potential energy, respectively. 
The deviation of the factor r\ from zero indicates, in principle, the overall quality of the 
variational wave function used [141 ] . In particular, for the wave functions used in our present 
calculations, the virial parameters n is in the range 1 ■ 10~ 18 — 5 • 10~ 19 , showing that our 
wave functions are highly accurate. 

Note that some of the bound state properties from Table II can be expressed as the linear 
combinations of other properties. For instance, by using the identity 

r 3 i = r 32 - r 2 i (10) 

one finds for the (r 31 • r 32 ) and (r 31 • r 32 ) expectation values 

(r3i-r 3 2) = ^((r3 2 2 ) + (r 3 2 1 )-(r 2 2 1 )) (11) 



and 



(r2i-r 3 2) = i((r 3 2 2 ) + (r 2 2 1 )-(r 3 2 1 )) (12) 



here and everywhere below the notation a • b stands for the scalar product of the a and b 
vectors, i.e. a-b = a x b x + a y b y + a z b z , where a = (a x , a y , a z ) and b = (b x , b y , b z ). Analogously, 
since (pi + p 2 + p 3 ) | \E r ) = 0, we have 

(Pi ■ Pi) = k(pl) + (P*) ~ (Pi)) (13) 
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where k) = (1, 2, 3). In the gradient form this equality takes the form 

(v, • v,) = + (vj> - (vl)) = -<4 V *> - <"^ v ?) + <-^ v '> ( 14 ) 

where (-|Vf) (i = 1, 2, 3) are the single-particle kinetic energies. In the general case, the 
both sides of the equalities Eqs.( !T3~]) - ffT4|) can be computed separately, i.e. such relations 
can be also used to control the overall quality of the variational wave function. The (pi • 
P2) expectation value can also be used as a measure of the dynamical electron-electron 
correlation in the two-electron ions. It should be noticed that there is an obvious difference 
between the (pi\l/ | p2^) and | Pi • P2 | expectation values. In fact, these two 
expectation values differ from each other by sign, i.e. (pi^/ | P2^) = — I Pi • P2 I 



A. Singular expectation values 

All expectation values mentioned above are regular, i.e. their analytical and/or numerical 
computation is relatively simple. A few expectation values from Table II, however, contain 
some singular parts. In such cases some additional explanations are needed. For instance, 
consider the expectation value (4~) = (4-)- In relative coordinates r 3 2,T 3 i,r2i the compu- 

r 32 r 31 

tation of this expectation value is reduced to the calculation of the following Hylleraas-type 
integrals 

T~2,i,i(a,b,c) = / / / exp(-ar 32 - br 31 - cr2i)r^ 2 r 31 r 2 idr 32 dr 31 dr2i (15) 

J JO J\r 3 2-r 3 i\ 

for many different sets of (a, b, c)— values. Formally, each of these integrals diverges, i.e. it 
does not exist as a finite expression. To make this integral finite one needs to introduce a 
small (and positive) cutoff parameter e 

/•oo roo |T32+r3i 

r- 2 ,i,iO> b,c;e)= / / exp(-ar 32 - br 31 - cr 2 i)r^ 2 r 31 r 2 idr 32 dr 3l dr 2 i (16) 

Je JO " , |r 3 2-r 3 i| 

Now, this integral is finite for e > 0, but diverges when e — > 0. For the finite integral 
r_ 2 ,i,i(a, b, c; e) defined by Eq.(IT6l) for e > we can write 

r_ a|1)1 (a,6,qe) = ^ (17) 

where 

r- 2 ,o,o(a, b,c;e)= / / exp(-ar 32 - br 31 - cr 2 i)r^ 2 dr 32 dr 31 dr 2 i (18) 

Je JO " , |r 3 2-r 3 i| 



The r_ 2 ,o,o(o, b, c; e) integral is represented as the sum of its regular (R) and singular (5) 
parts, i.e. r_ 2j0 ,o(a, b, c; e) = i?_ 2 ,o,o(«, b, c) + S- 2 ,o,o(a, b, c; e), where 

2 2 

<S-2,o,o(a, q e) = r~— f^( 1 ) ~ Zne ) = "F - v7s + ( 19 ) 
o + c v 7 o + c 

(o + c)ln(a + c) - (a + b)ln(a + b) 2 

#-2,o,o(M>c) = 2 — + 7^— ( 2 °) 

[tr — c z ) b + c 

= R -to,o( a > b > c ) + F -2,o,o(a, b, c) 



where ifi(n) is the digamma function 15] (or psi— function defined in Eq. (8.360) from 16| ) 
Note that -0(1) = — 7_e, where 7 £ ~ 0.577215 ... is the Euler's constant. The ^20 o( a , b, c 
term is the 'logarithmic' term which contains the ln(a + c) and ln(a + b) expressions, i.e. 

(a + c)ln(a + c) — (a + b)ln(a + b) , 



while 

2 



^-2,0,0(0, b, c) 



6 + c 

The second order derivative of the regular part is 

^ l(a ' 6 ' c) = + JbTcf (21) 

The expectation values of both sides of this equation are 

<* I R- W (a. b, c) I *> = <* I ^-1°^ !*> + <* I ^_ I *> (22) 

In the exponential basis the / 6H f c \a matrix elements correspond to the 5(r 32 ) delta-function 
multiplied by a factor of 47r. The corresponding expectation value is 

(* I #_ 2 ,i,i(a, 6, c) I ¥> = (* | ^ ~ J I *> + ^(S(r 32 )) (23) 

Analogously, the expectation value of the singular part is 

(* I 5_ 2 , 1 ,i(a,6,qe) | *> = -47r(5(r 32 ))( 7i; + Ine) (24) 

where (5(r 32 )) is the expectation value of the #(r 32 ) delta-function. 

As follows from the definition r_ 211 (a, b, c; e) = -R- 2 ,i,i(a, b, c) + S'_2,i,i(a, b, c; e) the sum 

r- 2 ,i,i (a, b, c; e) + S , _ 2 ,i,i (a, 6, c; e) = iL-2,1,1 (a, &, c) (25) 



is the regular expression and it has the finite limit when e — > 0. The limits of both sides of 
Eq.fl25D are 



V ) 

'32 



lim 

e-S-0 



(tf | 4- I *)e + 47r(5(r 32 ))(7s + Ine)] = 47r(5(r 32 )> + (* | i I ( 26 ) 



32 



32 



(i-> = 47r<5(r 32 )} + 



This expression is considered as the expectation value of the -3- operator, i.e. we have 

r 32 

- I *)r (27) 

32 '32 

Note that in contrast with the expectation values of regular operators the (4~) does not 

r 32 

coincide with the (\P I -3- I expectation value. The presence of the finite difference 

r 32 

47r(5(r 32 )) is typical for the expectation values of singular operators. For the (-5-) expecta- 

n 

tion value the analogous formula takes the form [17] 



'21 



lim 

e-S>0 



1 



21 



*} e + 47r(<5(r 21 ))(7 B + Ine) = 47r(5(r 21 )) + 



1 

^21 



*)r (28) 



Note also that in Eq. (l23|) the second order derivative is 



d 2 R-%fl{a, b,c) [(a + c)ln(a + c) - (a + b)ln(a + b)]bc 



dbdc 
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-4 



(c 2 - b 2 f 
[ln(a + b) + l}c [ln(a + c) + l]b 



(29) 



(b 2 - c 2 ) 2 (b 2 - c 2 ) 2 

It is easy to understand that this formula is computationally unstable, when c — > b. Indeed, 
both the numerator and denominator in this formula — > when c — > b. This produces some 
troubles in actual calculations. The transformation of this formula to the computationally 



stable form when c — > b can be found in 



17J. The most detailed analysis of various singular 



integrals arising in two-electron atomic problems can be found in 18 1. 



IV. CORRECTIONS TO THE TOTAL ENERGY 



The expectation values from Table II can be used to determine some actual properties 
of the helium atoms in the 2 1 S— state. Here by the 'actual properties' we mean some linear 
combinations of our expectation values which can be measured in modern experiments. 
The most important of such properties are various lowest order corrections to the total 
non-relativistic energies. Formally, by computing all possible lowest order corrections, e.g., 
relativistic and quantum electrodynamics corrections, mass corrections, etc, we must obtain 
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the exact agreement with the energies measured in high precision experiments. In this 
Section we consider a few such corrections to the non-relativistic atomic energies of the 
2 l S— state in the two-electron helium atom(s). 

First, let us evaluate the field component of the total isotope shift for the 2 1 5'— state in 
helium atoms. The field shift is related to the extended nuclear charge distribution which 
produces the non-Coulomb field at distances close to the nucleus. It is clear that the largest 
deviations between the Coulomb and actual potentials can be found close to the atomic 
nucleus, i.e. for distances r ~ r e < A < a 0l where r e = a 2 ao is the classical electron radius 
and A = crao is the Compton wave length. Here and below, a = 7.297352568 ■ 10~ 3 is the 
fine structure constant and ~ 5.29177249 • 10~ n m is the Bohr radius. Note that all 

3- 



numerical values for the physical constants used in this study were chosen from 



12|, 



The general th eory of the field shift has been discussed extensively in a number of works 



(see, e.g. 



19], 20[, [21 J and references therein). 



In our earlier work [22j we obtained the following expression for the field shift (in atomic 
units) of the bound S(L = 0)— states in light atoms and ions (with Q < 6) 

e{< = | QpMrHI^) = !«p,(0)fl 2 (iife = yO 4 ■ (SM)g)'t (30) 

5 e 

where Q is the nuclear charge and R is the nuclear radius. Our formula for the field shift E^ 
bllows from the well known expression obtained by Racah, Rosental and Breit (see, e.g., 
ig|). Its explicit derivation can be found in the Appendix. The parameter A and related 

factor £ = j in Eq. (l30|) describe the actual charge/proton distribution in the nucleus. 

5 

Also, in this expression r e = -^-^ = a 2 a ~ 2.81794093 fm (1 fm (fermi) = 1 ■ 10~ 13 cm) 
is the classical electron radius. In general, the nuclear radius R ~ r e and its actual value 
depends upon the total number of nucleons A in the nucleus (R ~ A^). In other words, the 
field shift formally corresponds to the a 4 — correction to the energy levels, i.e. to the second 
order relativistic correction. Let us evaluate the shifts for the 2 1 S— state in the 3 He 



and 4 He atoms. The nuclear sizes used in our computations were i?( 3 He) = 1.880 fm 23 1 



and i?( 4 He) = 1.6773 fm [24j (see also [25[ and |26|). We have also selected zero value for 
the parameter A, i.e. £ = 1 in Eq. (l3"0~|) . This means that the uniform (or r— independent) 
proton density distribution over the volume of the nucleus is assumed for each of the 3 He and 
4 He nuclei. This assumption produces the following values for the field shifts i?{|( 3 He) = 
1.66061-10 a.u. and £j^( 4 He) = 1.32120-10~ 8 a.u. Numerical recalculation of these values 
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to cm' 1 , MHz, eV and other units is straightforward. 

Another example is the lowest order (~ a 3 ) QED correction to the non-relativistic ener- 
gies of the 2 1 S— state of the helium atoms. For an arbitrary bound S(L = 0)— state in the 
°°He atom the closed analytical formula for the lowest order QED correction is written in 
the form (in atomic units) 



AE® = -Qa 3 



19 
30 



— 2 In a — In Kr 



<5(r 3 i)>+cr 



164 14 , 
1 In a 

15 3 



(31) 



where a is the fine structure constant, Q (= 2) is the nuclear charge and 5 is the total spin of 
two electrons. For the singlet states we always have 5 = 0. Also, in this formula In Kq is the 
Bethe logarithm 5j, which is represented in the form In Kq = In ko + 2 In Q, where In ko is the 
charge-reduced Bethe logarithm. Our current numerical evaluation of the charge-reduced 
Bethe logarithm for the 2 X 5— state of the °°He atom is In fc ~ 2.98011831. This value is not 
very accurate and must be improved in future calculations. Therefore, the value of Bethe 
logarithm for this state is In K rj 4.36641267. By using this value of In K and expectation 
values from Table II for the one finds that for the 00 He atom AE^ « 1.1890905 ■ 10~ 5 a.u. 
This is the lowest order (~ a 3 ) QED correction determined for the 2 1 S— state in the °°He 
atom. 

For the two-electron helium atoms with the finite nuclear masses we need to determine the 
corresponding finite mass correction (or recoil correction, for short). The recoil correction to 
the lowest order QED correction in the case of S(L = 0)— states in two-electron atoms/ions 
is represented in the form 27| (in atomic units): 



AE 



(3) 
M 



M 



-M + l M 



AE^ + 



4a 3 r 31 



3M 



+ 2-lna-41nfs: G 



14 



-a 



3 f 



(<5(r 31 )} 
1 , 1 



(32) 



M 3 MvrV 



r 31 



where AE^ is the expression from Eq. (l3ip . while M (M 3> 1) is the nuclear mass (expressed 
in the electron mass m e ). By using our expectation values from Table II and charge-reduced 
Bethe logarithm determined for the 2 1 S— state in the °°He atom one finds from Eq. (l32|) that 



A£g } ( 4 He) 



1.1891233 ■ 10~ 5 a.u. and AE$( 3 Ke) « 1.1891340 • 10" 5 a.u. To obtain a 



slightly better accuracy in the last equation we can use the expectation values from Table II 
determined for the two-electron He atoms with the finite nuclear masses. In these cases one 



11 



finds from Eq.(|32]) A^g^He) « 1.1891230 • 10~ 5 a.u. and AE^^Re) « 1.1891319 ■ 10~ 5 



a.u. 



For 



The last example which we want to consider here is related to the Vinti identity 
two-electron atoms/ions this identity takes the form 

(Pi " P2> = f (r 3 i " t 32 (±- + 4")> + ~<— ) = Q(^^> + ~<— > (33) 
2 rfi r§ 2 2 n 2 rgj 2 r i2 

The ( r3 ^^ 32 ) = | (^) + ((^) — (^)) expectation value can be computed either directly 
(see Table II), or with the use of the following relation 



<^>->-^> (34) 



since (pi ■ p 2 ) = \{v\) + \ {pi) ~~ \(p\) an d (Pi) = (pI) i n an y two-electron system. All 
expectation values mentioned in the right-hand side of this equality can also be found in 
Table II. Analogous relations can be derived for the ( r3 *a^ 21 ) expectation value. Note that all 
expectation values in the right-hand side of Eq.( )34|) are regular, while the expectation value 
in the left-hand side of Eq.f l34|) contains the difference of the two singular expectation values 
(^J 1 ) and (^J 1 ). As follows from Eq. (134jl such a difference of these two singular expectation 
values is a regular value. This means that the two singular parts cancel each other completely. 
In general, an accurate numerical coincidence of the ( rsi :f^ 2 ) expectation values computed 

r 31 

directly and with the use of Eq. (I34|) is another important test for our highly accurate wave 
functions. An analogous test can be used for the ( r31 3 r21 ) expectation value. It is also 

r 31 

interesting to note that for the singlet states in two-electron ions the C 31 / 32 ) expectation 
values slowly vary with the nuclear charge Q. Moreover, these expectations values have a 
finite limit (~ —0.1755) when Q — > 00. 



V. CONCLUSION 

We have performed highly accurate computations of the excited 2 1 S(L = 0)— states in 
the two-electron helium atoms: °°He, 4 He and 3 He. The total energies and a large number 
of bound state properties have been determined for the 2 1 S— state in the °°He, 4 He and 
3 He atoms to very high numerical accuracy. By using our highly accurate variational wave 
functions we have also evaluated the expectation values of some singular operators. To the 
best of our knowledge this work is the first extensive study of the bound state properties 
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of the excited 2 l S(L = 0)— states in two-electron helium atoms and most of the bound 
state properties computed in this work have never been evaluated in earlier studies. The 
knowledge of accurate expectation values given In Table II allows us to determine various 
lowest order relativistic, QED and mass corrections to the total energies of the 2 1 S— state(s) 
in the °°He, 4 He and 3 He atoms. 

In our extensive numerical calculations of the excited n}S— states of the He atom(s) we 
have found no additional complications for our highly accurate procedure which is based on 
careful optimization of many non-linear parameters at small and intermediate dimensions. 
Briefly, this means that 'instability problem' mentioned in the Introduction is not critically 
important for our method. As follows from computational results obtained in this study 
our methods can also be used for highly accurate computations of the excited states in two- 
electron atoms and ions. In fact, we have performed some of such calculations. Table III 
contains our preliminary results for the total energies of the 3 l S— and 4 1 S'— states in the 
°°He atom. As follows from Table III the total energies of the 3 1 S and 4 1 S' bound states 
converge very fast. Finally, our results from Table III for the the S 1 *?— and 4 1 S'— states 
in the °°He atom are the most accurate total energies ever obtained for these states. It 
shows a great potential of our method for highly accurate computations of bound states 
in two-electron systems. In our test calculations of the bound 5 l S— state in the °°He atom 
the short term booster function has not been constructed. Optimization of the three boxes 
for the non-linear parameters [8j allows us to produce a very compact (100-term) wave 
function which corresponds to the energy -2.0211768515145 a.u. The difference with the 
'exact' energy for this state is ~ 6.5 • 1CT 11 a.u. Very likely, that further optimization of the 
non-linear parameters in the wave functions with 200 basis functions, Eq.(J3]), will produce 
the variational energy which is very close to the known energy for this state. Analogous 
situation can be found for other highly excited S-, P— and D— states in the °°He atom. In 
future studies we want to develop the new optimization strategy which can be used to chose 
the non- linear parameters in trial wave functions for highly excited S— , P— and D— states. 
This will allow us to produce very compact and highly accurate wave functions for these 
states. Computations with the use of very large wave functions (> 3000 basis functions) can 
be avoided for highly excited bound states in two-electron atoms/ions. However, the current 
situation with highly accurate computations of highly excited singlet S— states (e.g., for the 
8 1 S—,9 1 S—, 10 1 S'— and higher singlet S— states) in the helium atoms and helium-like ions 
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is far from satisfaction for our method, since the Hylleraas variational expansion (see, e.g., 
4|) still provides better accuracy for such states (if comparable numbers of basis functions 
are used in both methods). 

Note that our methods allow one to determine the total energies and many other bound 
state properties to very high numerical accuracy which is quite comparable and even bet- 
ter than analogous accuracy achieved for the ground states. The overall convergence rates 
observed in calculations of many bound state properties of the excited 2 1 S(L = 0)— states 
in two-electron helium atoms are relatively high. Briefly, this means that our optimiza- 
tion strategy used to optimize the non-linear parameters in highly accurate wave functions 
works very well for both ground and excited states in two-electron atoms/ions. This also 
contains the answer to the third problem mentioned in the Introduction that is the conver- 
gence rate for most of the bound state properties of the excited states is not substantially 
lower than for analogous properties of the ground state. Furthermore, the same optimiza- 
tion procedure has been applied to highly accurate computations of the 2 1 P(L = 1)— and 
2 3 P(L = 1)— states in the He atom(s). The explicit form of the trial variational wave func- 
tion for the states with L > 1 is given in [8|. The best variational (total) energies obtained 
in our computations for these states in the °°He atom are: -2.1238430864981013590742 a.u. 
and -2.13316419077928320510251 a.u., respectively 29[. Numerical uncertainties in these 
values can be evaluated as ~ 3.5 -10 -20 a.u. In other words, these energies and corresponding 
wave functions are significantly more accurate than known from the modern literature (see, 
e.g., b|). 

Our highly accurate method can also be applied to highly accurate calculations of the 
weakly-bound singlet and triplet excited states in two-electron atoms and ions. These are 
the n 1,3 S—,n 1,3 P—,n l,3 D— states with n > 2. Results of these and other similar computa- 
tions indicate clearly that clusterization of the wave functions plays a very important role 
for the excited states in two-electron atomic systems. The effect of clusterization for the ex- 
)onential variational expansion was originally discovered in 30] . For our two-stage method 
8] the clusterization of variational expansion, Eq.([3]), means the presence of two following 
things. First, the short term (booster) wave function produces almost 'exact' energy al- 
ready for N = 200 - 300. For instance, for the 3 1 S'(L = 0)— state the booster function 
with only iVo = 200 exponential basis functions allows one to obtain the total energy E = 
-2.061271989738157 a.u. (current value), while the 'exact' energy of this state in the °°He 
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atom is -2.061271989740911 a.u. Second, optimization performed at the second stage of our 
method produces a number (usually three) of 'optimal' parallelotops which are used later to 
chose parameters ctj,/^ and 7$ (where Nq < % < N) in Eq.Q. In actual calculations of the 
excited n 1 ' 3 ^— and n 1,3 P— states in the He atom(s) for n > 3, such an optimization of three 
parallelotops can be performed for TV = 600 in the total wave function. It produces almost 
the exact value of the bound state energy. Moreover, the optimal parallelotops for these 
bound states have almost degenerated structure. Briefly, this means that such parallelotops 
generate the wave function, Eq.([3]), in which all parameters 7, either equal zero exactly, 
or very close to zero. But this means that the role of electron-electron correlations (or r 2 i 
coordinate) rapidly decreases for highly excited bound states in two-electron ions [30]. In 
respect with this, the difference between the total energies of the triplet and singlet states 
with large angular moment L in He-like atoms and ions rapidly (almost exponentially) de- 
creases, if the value of L increases. This remarkable property of the exponential variational 
expansion, Eq.Q, in applications to Rydberg states has been observed in all helium-like 
atoms and ions. It can be used to simplify future highly accurate computations of highly 
excited bound states in two-electron systems. 
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Appendix 

The formula obtained by Racah, Rosental and Breit for the field shift is 

fs Anal r/ \ 7+1 m s (2QR\^ SR 

M = ~Q~ ' ( TeN ^ ' f(27+T) ' ^ ' (~a^~J ll (35) 

where Q is the nuclear charge, R is the nuclear radius and 7 = a/1 — a' 2 Q' 2 . For the helium 
atoms we have 7 ~ 0.9998935. Therefore, to a very good accuracy one can assume that 
7 f=s 1. The factor B(j) in Eq. (l35|) is 

Bh) = 7 £ r w - (36) 

KU (2 7 + l)(2 7 + 3) 5 V ; 



It corresponds to the uniform distribution of the proton density over the volume of the 

SR 
R 
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nucleus. The ratio §§■ in the formula Eq. (1351) equals unity, if the field shift is determined in 



respect to the atom with a point nucleus. Now, in atomic units the formula Eq. (l35j) takes 
the form 

where r e = -^-^ = a 2 a is the classical electron radius. Finally, we obtain the following 
formula for 

E J M = —a ■ 5{r eN ) ■ f-j (38) 

This formula can slightly be modified to include other possible distributions of the proton 
density in the nuclei. The formula given in the main text contains an additional factor 
fW = ^ e f° rmma Eq. (j38p with such a factor (i.e. the formula Eq. ( I30|) from the 

5 

main text) also provides the correct answer in those cases when the uniform distribution 
of the proton density over the surface of the nucleus is assumed (in this case A = 1). 
Furthermore, variations of the parameter A in Eq. (l30|) (0 < A < 1) allow one to describe 
other possible distributions of the proton density in light atomic nuclei with Q < 6. 
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TABLE I: The total non-relativistic energies E and (r2i),(#3i), v^i expectation values computed 
for the excited 2 1 S(L = 0)— states of the °°He atom (in atomic units). 



N 


£(°°He) 


( r 2l) 




^31 


1500 


-2.145 974 046 054 416 890 40 


5.26969620234182140 


1.3094607823 


-2.0000000848 


2000 


-2.145 974 046 054 417 174 01 


5.26969620234181075 


1.3094607799 


-1.9999999787 


2500 


-2.145 974 046 054 417 298 64 


5.26969620234180607 1.3094607811 


-2.0000000487 


3000 


-2.145 974 046 054 417 342 16 


5.26969620234180443 


1.3094607802 


-1.9999999798 


3500 


-2.145 974 046 054 417 372 44 


5.26969620234180331 


1.3094607801 


-1.9999999784 


4000 


-2.145 974 046 054 417 385 25 


5.26969620234180281 


1.3094607802 


-1.9999999977 


4200 


-2.145 974 046 054 417 391 41 


5.26969620234180258 


1.3094607804 


-2.0000000089 


oo( a ) 


-2.145 974 046 054 417 415(10) 






-2.0 




-2.145 974 046 054 419(6)^ [4] 










-2.145 974 04(1)W [3j 









(a) The asymptotic value of the total energy (in a.u.) 
( 6 )The best variational results known from earlier calculations [4] and [3] in atomic units. 
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TABLE II: The bound state properties X computed for the excited 2 1 S(L = 0)— state in the °°He, 
4 He and 3 He atoms (in atomic units). 



X 


°°Hc 


"He 


3 He 




-2.14597404605441739141 


-2 14567858758314906950 


-2.14558192369821233185 


(T) 


2.14597404605441739081 


2 14567858758314906890 


2.14558192369821233125 


(V) 


-4.29194809320883477957 


-4.29135717516629813451 


-4.29116384739642465859 




4.57M0" 19 


7.8674- 10" 19 


9.1373-10~ lfl 


(''si 2 } 


0.1437248133044633 


0.1436882345740985 


0.1436762703243981 


<>'2i 2 } 


4.1469390197898792 


4.1458017822058179 


4.1454297481103379 


i r 3l) 


1.13540768612560041 


1.13525123516046493 


1.13520005025470815 


i r 2~l) 


0.24968265239356710 


0.24964776547556185 


0.24963635362241032 


< 1 > 


0.39836585186133074 


0.39819906128376390 


0.39814450960354474 


^ T21J-31 ) 


0.34063384586003005 


0.34053958784463326 


0.34050875587573647 




0.56186146745960207 


0.56170295659588612 


0.56165110690333470 


z L r 3i r '.n 


-0.625372574195182710 


-0.625157759769931583 


-0.625087487579971712 


/ 121 *22 \ 

r 31 


-0.0576687311323827964 


-0.0575321421896994047 


-0.0574874624526176361 


/r;ii-r 2 i \ 


1.19307641725798262375 


1.19278337735016376067 


1.19268751270732578757 


(r™ ) 
V 31/ 


2.973061 134389489243 


2.973491187710339638 


2.973631893971735653 


V 21/ 


5.269696202341802818 


5 27045091 723581 4928 


5 270697843066571 Q09 


/ r 2 \ 
v 31/ 


16.08923324404050295 


16.09391391074115758 


16.09544547281748570 


V 21/ 


32.30238037187441050 


32.31159744722066081 


32.31461335707651003 


(rh ) 
V 31/ 


108.06042106146845 


108.10756711191265 


108. 12299522961221 


V 21/ 


224.03043954850232 


224. 12604278314834 


224. 15732795934735 


V 31/ 


825.75317813369571 


826.23293858862828 


826.38995080810075 


(rt) 

V 21/ 


1737 41 5292784614 


1 738 402036368951 


1 738 7949681 94081 




0.557144578327034669 


0.557137512037319818 


0.557135200602530583 


T"21 


-0 01 46570433571 77936 


-0 01464331 5290865667 


-0 01 46388941 1 0935795 


W / 


2.49080283242230255- 10~ 2 


9 490 79 9 7 1 9 9 1 3 1 9 9 4- 1 0~ 2 


9 4907894973531 3599' 10~ 2 

^.-t!y\J 1 (1.1 1_ f T }'.) ' } 1 r > '.) -J o ±\J 


\ 2 V 1/ 


1 0799870230972085694 


1 0796915861654766111 


1 0795949387331539546 


\ 2 V 3/ 


2.1554779104737081945 


2.1548473559242290248 


2.1546410836274219596 


(Pi 1 Pa) 


-9 1 5547791 04737081 945 


-2.1548473559242290248 


-2.1546410836274219596 


(pi • P2) 


Q c ifl^Hfi44 l Q9Q1 0^(SXQ4- 1 n -3 

,/.:JUv>OU 1 1 L.I £..) L\.!:>: IllOyi 1U 


9 4641 835939758096Q7310~* :i 


q 4^1 90fi i tsi 1 1 4n c in^^ c in.i n -3 


( r 21 ' r 3l) 


16.1511901859372052475 


16.1557987236103304049 


16.1573066785382550123 


{^*31 * ^32) 


-6. 195694189670229434- 10~ 2 


-6.188481286917282672-10 -2 


-6.186120572076931285-10~ 2 


\°V l 31)/ 


1 '^094607809 


1.3089223177 


1.3087461812 




-1.999999997659 


-1.999725848502 


-1.999636146310 




-2.0 


-1.999725850875267686059 


-1.999636157582479619549 


<«(r 2 i)> 


8.648433654- 10~ 3 


8.64511990M0- 3 


8.644035949- 10~ 3 


V21 


0.4999998171 


0.4999994703 


5.0000004784 


v 2\ 


0.5 


0.5 


0.5 


(fel) 


0.1755555002 


0.1753991714 


0.1753446276 


('2l 3 >« 


-0.04073283492915455 


-0.04069706227378580 


-0.04068536285394055 


(r^h 


-39.23799206891561 


-39.21960245645624 


-39.21358724048685 


(r2l) 


0.06794658607018129 


0.06794071686462781 


0.06793879648060461 


<>'31 3 > 


-22.7828226037360592 


-22.7711995110946440 


-22.6896828517925066 



(a) The exact value. 
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TABLE III: The total non-relativistic energies E (preliminary results) of some excited n 
0)— states (n = 3, 4) of the °°He atom (in atomic units). 



N 


E(3 1 S -state) 


S(4 1 5-state) 


1500 
2000 
2500 
3000 
3500 
3800 


-2.061 271 989 740 906 234 
-2.061 271 989 740 907 482 
-2.061 271 989 740 907 999 
-2.061 271 989 740 908 266 
-2.061 271 989 740 908 387 
-2.061 271 989 740 908 430 


-2.033 586 717 030 718 055 
-2.033 586 717 030 722 441 
-2.033 586 717 030 723 921 
-2.033 586 717 030 724 515 
-2.033 586 717 030 724 857 
-2.033 586 717 030 725 074 


ooW 


-2.061 271 989 740 908 48(5) 


-2.033 586 717 030 725 20(7) 


f 4 ] (") 


-2.061 271 989 740 911(5) 


-2.033 586 717 030 72(1) 



( a )The asymptotic value of the total energy (in a.u.) 
( 6 )The best variational results known from earlier calculations [4] in atomic units. 
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